Universal hadronization condition from RHIC to LHC 
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We obtain a detailed description of all available hadron multiplicity yields in geometrically central 
Pb-Pb collisions at LHC measured in the rapidity interval \y\ < 0.5. We find that the hadronization 
of the fireball at LHC occurs at nearly identical intensive physical bulk conditions for all geometrical 
centralities similar to those seen already at RHIC. 
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Introduction and motivation: We extend the suc- 
cessful description of central rapidity particle yields in a 
single freeze-out model [H Q to characterize the physical 
properties of the hadronizing fireball. We consider, as an 
example, a supercooled quark-gluon plasma (QGP) dis- 
integrating into hadrons which can scatter but preserve 
the final stable particle abundance. Therefore hadron 
particle multiplicities directly characterize the properties 
of the fireball. Final state hadrons are thus produced 
according to the accessible phase space with otherwise 
equal reaction strength. Accordingly, the particle yields 
are described by the chemical non-equilibrium statistical 
hadronization model (SHM) 

In this SHM implementation within the SHAREv2.2 
program, the yields of particles are given the chemi- 
cal freeze-out temperature T and overall normalization 
dV/dy (paralleling the experimental data available as 
dN/dy). We include phase space occupancies 7g,7s for 
light {q = u,d) and strange s quark flavors, respec- 
tively, and we account for the small asymmetry between 
particles and anti-particles by fugacity factors Xq, Xg 
and the light quark asymmetry A/3. These parame- 
ters enter the distribution function as fi(s,T,ji, Xi) = 
l/(7riA±ie^/^ + S) for flavor i, where S = -1,0, +1 
for bosons, Boltzmann distribution and fcrmions respec- 
tively. Wc refer the reader to Section 2 of Q for further 
discussion of the above parameters. 

Our discussion addresses, evaluates and compares the 
LHC Pb-Pb experimental results available at ^/sWn — 
2760 GcV and RHIC Au-Au at y/s^N = 62.4 GeV. Wc 
show that the chemical non-equilibrium SHM works at 
LHC in the 0-20% geometric centrality. This is so, since 
in the chemical non-equilibrium SHM approach, we allow 
quark pair yield parameter 7^ > 1 for light quarks; this is 
the key difference from the simpler equilibrium SHM. The 
rationale for 7q 7^ 1 originates in the high entropy density 
of QGP at hadronization compared to the hadron phase 
space in which the color degree of freedom is frozen. In 
case that the fireball disintegrates faster than the time 
necessary to equilibrate the yield of light quarks bound 
in hadrons, the value 7^ > 1 must arise. 

One can consider a dynamical volume growth as a path 
for absorbing the excess entropy of QGP source. This 
requires a much longer lifespan of the particle source 



than is supported by HBT data [SL l6|. The observed 
LHC total lifespan (r/ ~ lOfm/c [6|]) favors very fast, 
or sudden, hadronization 0, 01 ■ In this situation, chemi- 
cal non-equilibrium approach must be applied also to the 
light quark abundance, introducing the light quark phase 
space occupancy 7n- This proposal made for the high en- 
ergy SPS data [9|, llO| helped also to improve the under- 
standing of RHIC200 hadron rapidity yield results [llj 
and allowed a consistent interpretation of such data 
across the full energy range at SPS and RHIC200 12|. 

SHM fits, that arbitrarily set 7, = 1 H, 14|, (equilib- 
rium SHM) , describe hadron yields at LHC with relatively 
large total x^- This chemical equilibrium SHM disagrees 
at LHC across many particle yields, but the greatest is- 
sue is the 'proton anomaly', which makes it impossible 
to fit the p/tt = 0.046 ± 0.003 ratio ^ along with the 
multi-strange baryons S and fl. 

We will show, that common intensive QGP bulk prop- 
erties arise, nearly exactly equal to those found at RHIC, 
for all four geometric collision centralities we analyze. 
We will discuss in depth the main extensive bulk prop- 
erty difference we find, that is the entropy dS/dy growth 
with energy and centrality. 

Fit to most central collisions: Within the chemi- 
cal non-equilibrium SHM wc have, allowing for baryon- 
antibaryon asymmetry, seven independent statistical 
model parameters reduced by two constraints: a) The 
isospin fugacity factor A/3 is constrained by imposing the 
charge per baryon ratio, ((Q) - (Q))/((B) - (B)) ~ 0.38, 
present in the initial nuclear matter state at initial in- 
stant of the collision; b) For each value of Xq, strangeness 
fugacity As is evaluated by imposing the strangeness con- 
servation requirement (s) — (s) ~ 0. Considering the 
particle-anti-particle symmetry at LHC, the four key pa- 
rameters are the hadronization volume dV/dy, tempera- 
ture T and the two phase space occupancies 7^ and 7s. 
The 7th parameter is light quark fugacity Xq. 

We use as input to our fit the hadron yield data in 
0-20% centrality bin as presented in Ref. |l5|, where a 



fit to this data set for the case of chemical equilibrium 
(7s = 75 = 1) is shown. For comparison and demonstra- 
tion of method compatibility, the chemical equilibrium 
model fit (dashed lines in figured]) is shown, with a large 
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Data: ALICE, 0-20% (preliminary) 
Model: Non-equilibrium SHM 
Model: Equilibrium SHM 
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FIG. 2: (color online) Experimental data (full symbols) and 
model predicted (open circles) for particle ratios as a function 
of centrality. See text for discussion of data and results. 



FIG. 1: (color online) The non-equilibrium SHM fit is indi- 
cated by (blue) solid horizontal lines overlaying for all the 
LHC-ALICE (preliminary) data available in 0-20% central- 
ity bin (red squares). Chemical equilibrium fit is indicated 
by (green) dashed lines with model parameters presented in 
parentheses. 

= 34.8. In both approaches, we fit the same data, 
the decrease in by factor nearly 5 is due to chemical 
non-equilibrium i.e., ^ li 7,s 7^ 1- We determine the 
best hght quark fugacity factor A, — 1.00359, which cor- 
responds to baryo-chemical potential fis — 1.5 MeV, and 
we apply strangeness and charge per baryon conservation 
by fitting them as two additional data points. The result 
of our 0-20% centrality bin fit is shown in figured] and in 
top section of the third column of table [11 W e compare 
our present results to our recent analysis [ig of Au-Au 
coUisions at ^/snn = 62.4 GeV at RHIC62, shown in the 
2nd column in table HI 

More peripheral centralities: Wc extend our study 
to more peripheral collisions at LHC using much smaller 
data set, complemented by two assumptions as follows: 

a) We consider the three ratios K*°/K^, A/tt = 
2A/(7r~ -j-TT"'') and as presented in [l3| and which 
we show in figure [5] with full symbols. We fit these ratios 
in three centrality bins, 20-40%, 40-60% and 60-80%, in 
which K*/K and A/tt have experimental data point. We 
take average of two neighboring (/i/K data points in order 
to use this ratio as input to our fit in the intermediate 
centrality. This is consistent with the claim that 0/K is 
constant over all centralities and has been claimed inde- 
pendent of centrality in [l7| . 

b) To obtain overall normalization, we complement the 
ratios with charged particle rapidity density dNct/dy. 
Based on our fit of 0-20% centrality data, we see that the 
ratio of charged particle rapidity to pseudo-rapidity den- 
sity is {dNcti/ dy) / {dNchI drf) = 1.115. We multiply data 



from [3l by this factor, and use the resulting dNch/dy 
as an additional data point, that determines the value 
of fireball volume dV/dy. The input multiplicity data is 
presented in the figure [S] 

c) We include the yields of tt^, and p^, as presented 
in [l^ and shown in figure |31 Similarly to (/)/K ratio, we 
use as input averages of yields in two neighboring cen- 
trality bins. 

Wc fit the three centrality bins of LHC2760 data us- 
ing the eight data points and the two conservation laws, 
which fix As and A/3. We thus have five degrees of free- 
dom as is seen in table U in 4th, 5th and 6th column. All 
studied centralities show reasonable x^/ndf < 1.1. We 
find for all four LHC centrality bins remarkably similar 
statistical parameters. 

We compare the outcome of the fit showing the three 
input ratios K*°/K~, A/tt and 0/K~, and p/tt in figureO 
All data are well fitted including p/tt, which we evaluate 
from the individual fitted yields of p and tt from [l^ . 

The predicted particle yields normalized by iVpart/2 
are shown in figure |31 where considering particle-anti- 
particlc symmetry only one of the isospin multiplets is 
shown. The 7r,K and p yields are fitted values, whereas 
the other particle yields are predictions. We indicate in 
figure [3] on the right edge the experimental input for the 
0-20% bin with an offset to assure visibility of the small 
differences between fit and experimental data. We show 
the input multiplicity {dNch/dy)/{Npai-t/'2) and fit result 
in figure [3l Both are overlapping exactly for the three 
peripheral bins since this is the most precise input data. 

Hadronization conditions: Despite a change by a fac- 
tor of 45 in reaction energy comparing RHIC and LHC, 
and the wide range of centrality, the only quantity among 
statistical parameters shown in table H] that significantly 
changes is dV/dy. This suggests, that we should look 
closer at the intensive bulk physical properties of the fire- 
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TABLE I: Top section shows chemical non-equilibrium SHM fit parameters dV/dy, T, 7^, 7^ and Xtotai with ndf (number data 
less number of parameters) obtained in each centrality bin. Errors are a fit-stability estimate obtained with yield shifted 
within experimental error, the underlying statistical fit error is negligible. Bottom section presents fireball bulk properties 
in each bin: energy density e, pressure P, entropy density a, strangeness per entropy content s/S, and entropy at LHC2760 
compared to RHIC62, Slhc/Srhic. The centrality defining number of participants Npart values are adopted from [20| . 





RHIC62 


LHC2760 


Centrality 


0-5% 


0-20% 


20-40% 


40-60% 


60-80% 


(iVpart) 


346 


308 


157 


68.8 


22.6 


dV/dy [fm=*] 


853 


2455 ± 146 


1169 ±9 


406 ±3 


102 ± 7 


T [MeV] 


139.5 


138.6 ± 1.1 


137.6 ± 0.03 


140.5 ±0.04 


143.2 ± 0.08 


7? 


1.58 


1.627 ±0.007 


1.633 ± 0.0002 


1.616 ±0.003 


1.60 ±0.02 


7s 


2.24 


2.04 ±0.04 


2.01 ±0.12 


1.83 ± 0.08 


1.70 ±0.09 


Xtotai/ndf 


0.38/5 


7.40/8 


2.93/5 


3.58/5 


5.43/5 


e [GcV/hV] 


0.493 


0.466 ±0.018 


0.441 ±0.012 


0.488 ±0.010 


0.536 ± 0.025 


P[MeV/fm^] 


82.0 


79.1 ± 2.8 


75.5 ± 1.5 


82.2 ± 1.3 


90.2 ±4.0 


a [im-^] 


3.40 


3.23 ±0.11 


3.07 ±0.07 


3.36 ± 0.06 


3.65 ±0.13 


s/S 


0.0322 


0.0296 ± 0.0002 


0.0289 ± 0.0014 


0.0277 ± 0.0009 


0.0267 ±0.0011 


•S'lhc/5'rhic 




3.05 


2.66 


2.18 


1.52 
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FIG. 3: (color online) Predicted particle yields per participant 
pair as a function of centrality. Open symbols represent our 
model predictions; lines guide the eye. We show experimental 
data for the 0-20% bin as full symbols, with an offset. 



ball: The emitted hadrons not only carry away from the 
fireball the above discussed charge, baryon number or 



strangeness, but also e.g. the thermal energy, dE/dy ob- 
tained by summing the energy content of all produced 
particles, observed and predicted. 

The bulk thermal energy density at hadronization e 
defined by e = {dE / dy) / {dV / dy) is of direct interest. 
Similarly, we evaluate the entropy dS/dy, pressure P and 
total yield of strangeness ds/dy = d{qs ± qs)/2dy. These 
properties of the fireball at hadronization are shown in 
the bottom section of the table HI where for simplicity we 
omit in the first column the symbol d/dy. 

As one can see by comparing the second and third 
columns of table HI the intensive properties of the RHIC62 
fireball at hadronization i.e. s,P,a and s/S, are practi- 
cally identical to the here evaluated case LHC2760, and 
this continues across all considered centralities as is seen 
in the 4th, 5th and 6th column of tableH) We show a com- 
parison of e, P as a function of centrality between LHC 
(solid symbols) and RHIC (open symbols) in the bottom 
part of figure The difference between LHC and RHIC 
can be easily attributed to the fit uncertainties, since the 
intensive quantities are proportional to a high power of 
statistical parameters. 

Hadronization volume dV/dy does not characterize 
the early stage of a collision, this information is avail- 
able in the entropy content at hadronization dS/dy cx 
dV/dy which presents a more accurate view of the pre- 
hadronization processes that created the fireball. For an 
ideally flowing and expanding QGP, most of the observed 
entropy yield dS/dy of a fireball is created early in the 
collision. In the top panel of figure |4] we sec that at 
LHC, dS/dy oc N^J^j''^, the rise is faster than linear. For 
comparison note that at RHIC, the entropy yield rises 
almost linearly with Npart- Last row of table H] shows 
the enhancement of entropy at LHC compared to RHIC, 
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FIG. 4: (color online) Top panel shows entropy content of the 
fireball at LHC (full symbols) and RHIC (open symbols) as 
a function of centrality. Bottom part shows pressure P and 
energy density e at hadronization with the same symbols for 
LHC and RHIC as in the top panel, for values see table HI 



S'lhc /'S'rhic • The enhancement decreases as a function 
of centrality from ^ 3 to ~ 1.5, which implies additional 
entropy production mechanism proportional to centrality 
at LHC. 

Strangeness per entropy s/S = {ds / dy) / {dS / dy) is of 
particular interest in the source fireball since both en- 
tropy and strangeness yields are nearly preserved in the 
hadronization process, but the production of strangeness 
occurs after most of the entropy is created. Up to 
a well studied proportionality factor, s/S is the ratio 
of strange quark abundance to total quark and gluon 
abundance which is making up the entropy in the bulk. 
Therefore, s/S measures the degree of chemical equi- 
libration attained in the QGP We observe a constant 
value of s/S ~ 0.03 (see table H]), which is in agreement 
with theoretical expectations for the strange quark mass 
nis ~ 100 McV 

Comments and Conclusions: The LHC2760 exper- 
imental environment has opened a new opportunity to 
investigate in detail the hadron production mechanisms. 
Precise particle tracking near to interaction vertex in the 
ALICE removes the need for off-line corrections of weak 
interaction decays, and at the same time, vertex tracking 
enhances the efficiency of track identification, increas- 
ing considerably the precision of particle yield measure- 
ment 
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15l |. All LHC experimental results used in the 
present work were obtained in this way by the ALICE 
experiment for Pb-Pb collisions at ^/snn = 2.76 TeV, 
limited to the central unit of rapidity interval —0.5 < 



y < 0.5. 

In this new experimental environment we show the 
necessity to introduce the final state hadron chemical 
non-equilibrium, which describes well all experimental 
results obtained in the Pb-Pb collisions at y^s]v7V = 
2.76 TeV from LHC. As figure [1] shows, jq ~ 1.6 (non- 
equilibrium of light quarks) allows to describe the ratio 
p/tt = 0.046±0.003 together with yields of mufii- 

strange baryons S and fl. 

Another approach to describe the data including the 
'anomalous' proton yield at LHC involves chemical equi- 
librium hadronization at relatively high T followed by 
hadron interactions 



23|. We note the chemical equi- 
librium SHM yields at hadronization in figure [U which 
overpredict proton yield, and at the same time underpre- 
dict both S and fl. Any alternate data explanation must 
come to terms with this situation, thus it must deplete 
protons and enhance both S and Q, and at the same time 
the ratio p/tt must remain practically constant . This is 
difficult, as we now discuss, looking closer at the results 
ofRcf. [13, [H: 

• we see in figure 1 of (2^ that if and when equilibrium 
style hadronization occurs and leads to high T the neces- 
sary post-hadronization reactions deplete protons and S 
and enhance of Q. This means that the already too small 
a yield of S is further depleted and disagrees gravely with 
experiment. 

• the measured p/tt ratio in the 0-5% centrality bin can 
be made consistent with post-hadronization proton-anti- 
proton annihilation [2^. This fine tunes model parame- 
ters and as a result for the 20-30% centrality bin Ref. [i^l 
reports increased p/tt ~ 0.058. 

• the model predicts for peripheral collisions yet less an- 
nihilation and thus a p/tt ratio approaching equilibrium 
SHM value, which is twice as large as experiment. While 
experiment for p/tt seen in figure [2] is a constant for all 
centralities, Ref. 23| thus predicts a rapid variation by 
about factor of two. 

These arguments lead to the conclusion that post- 
hadronization interactions are inconsistent with the ex- 
perimental data of baryon yields at LHC. On the other 
hand, our chemical non-equilibrium SHM at LHC pro- 
duces a high confidence level fit x^/ndf = 7.4/8 < 1. 
Prior SPS and RHIC data analysis [Til, [13, [H has al- 
ready strongly favored chemical non-equilibrium variant 
of SHM. The implied sudden hadronization picture is per- 
fectly consistent with the anisotropic flow of quarks lead- 
ing to the final hadron momentum distribution azimuthal 
asymmetry (see e.g. [ist). 

Moreover, we find that LHC and RHIC results are 
quite consistent in our approach, we obtain the same 
hadronization condition (e, P, a) at LHC as previously 
reported at RHIC, which in turn agrees with high en- 
ergy SPS [I4I. The energy density of hadronizing mat- 
ter is 0.50 ± 0.05 GeV/fm^, which is about 3.3 times 
the energy density of nuclear matter, the pressure is 
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P = 82 ± SMeV/W = (158 ± AMeV)^, as is seen in 
the bottom part of figure [H and which has been pro- 
posed in [2^. The bottom part of table U also shows that 
the entropy density is constant: a — 3.35 ± 0.30 fm~'^ for 
both experiments and all centralities. 

These typical QGP properties, including s/S ^ 0.03, 
mean that at LHC, the source of hadrons is a chemically 
equilibrated strangeness saturated QGP fireball. Further- 
more, the universal hadronization condition cannot be 
viewed anymore as being due to successive particle emis- 
sion, or to proceed via equilibrated hadron gas phase. 
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